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Abstract: We discuss several issues related to recent work on high energy evolution. In particular 
we show that the Hilbert space for action of the operator of the evolution can be conveniently 
defined by the space of weight functional introduced recently in |28| ]. These weight functionals 
automatically implement the correct properties of the charge density correlators, thus eliminating 
the need in explicit introduction of the Wess Zumino term. We also discuss various aspects of 
Dense Dilute Duality in the toy dipole model. 
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1. Introduction 



The ideas of gluon saturation]!]] have been vigorously pursued in the last ten years or so [§, |3|, ||[ ||, 
H|, ^, p. The result of these efforts is the set of evolution equations which describes the change of 
the scattering amplitude at high energy in the regime where the hadronic wave function resembles 
a dense gluonic medium - the so called JIMWLK equation || ^. [?], ^] (for reviews see ||, 1C, 11, |12||). 

The latest developments in this dynamic field are related to what is sometimes called Pomeron 
loops, or Fluctuations or alternatively Saturation effects in the projectile wave function [13, 14, 
1§ [fi| 0, 0, 0, 0, |2|, |2| |24], 



2§ [27], |28|, |28|, |30|, |31|, |32 

to these effects clearly must be important for the preasymptotic regime, but should also affect the 
more exclusive observables at very high energy. The last year or so has seen a lot of activity in 
this direction. Much has been understood during the last year, but many questions remain open. 



35]. The corrections due 



In particular the structure of the "Hilbert space" on which the evolution kernel acts [24| as well 
as the precise physical meaning of the Dense-Dilute Duality (DDD) introduced in [p6| l are being 
actively debated (p8L 31, 32, 



34] 



This paper presents several remarks on these issues. The remarks are not necessarily causally 
or logically connected, but they all pertain to the questions mentioned above. 

As regards the first issue, one of the questions is the following. The evolution is represented 
usually as acting on the weight function W[p(x, x~)], where the necessity of introducing the longi- 
tudinal coordinate x~ stems from the desire to represent quantum noncommuting operators p a (x) 



in terms of commuting " classical fields" . It was shown in |24|] that this is achievable by endowing 
the quantum operator with an additional "ordering" variable. One also needs the weight functional 
W to ensure that the correlators of the type (p(x^)...p(x n )^) = J D[p]W[p]p(x^)...p(x~) depend 
only on the ordering of the variables xj rather than on their values, and that an exchange of two 
such variables must lead to the appearance of the SU(N) structure constant. Namely if in the 
correlation function the two values of the ordering variable x^,xjare adjacent to each other and 
x7 L > xj , the following property must hold 

D[p]W[p]...p a (x-)p\xj)... = J D[p}W[p]...p a {xj)p b (xr)... + if abc J D[p]W[p]...p c {xT)... 

(1.1) 



It was also shown in [24] that these requirements can be satisfied by a "weight function" W 
which contains a Wess-Zumino term. The explicit representation of such a Wess - Zumino term is 
awkward, since it has to be local in transverse space and has to allow for an arbitrary representation 
of the colour charge at every transverse position p4 



On the other hand in [25] we have introduced a different representation of the weight functional 



where 



R ab (x) 



W[p a ] = E[R]5[p(x,x-)} 



Pexp{ / dx'T 



;} 



ab 



(1.2) 



(1.3) 



5p c (x, x 

and £ - an arbitrary functional. A similar representation was introduced earlier in the framework 



of the dipole model[13], with the crucial difference that the exponential in functional derivatives 



was truncated at a low order in the Taylor expansion. Although the formula eq.(1.2) has a very 
natural interpretation from the point of view of partonic content of the wave function, the question 
whether it satisfies the property eq. (|l.l|) has not been addressed before. Our observation is that 
the Hilbert space for the action of the high energy evolution kernel is indeed naturally given by 



the vector space eq.( |1.2[ ). In particular we show that the necessary condition eq.(l.l) is satisfied 
by any element of the space defined by eq. ( |1.2[ ) . We thus note that this form constitutes a simple 
alternative to the Wess-Zumino term as far as the correct representation of the charge density 
commutation relations is concerned. This discussion along with some additional remarks which 
are mostly meant to clarify and explain in more detail our earlier results ]^4|, 25|, 26, ^] is contained 
in Section 3. 

In Section 4 we discuss a dipole toy model with no transverse coordinates. We ask the question 
what is the DDD transformation and what is the analogue of the Seld-Duality condition in terms 
of the dipole degrees of freedom. By applying the DDD transformation we construct the dual 
to the Mueller's dipole evolution. In the framework of this model we explicitly construct a self 
dual evolution which is valid both in the dilute and dense limits. In the intermediate regime the 
model is incomplete and misses correct treatment of multiple dipole emissions. This is reminiscent 
of the self dual variant of the high energy evolution suggested in [31] and discussed recently in 
J34j[ . The value of this discussion is mainly in the simplicity of the model and transparency of the 



mathematics involved which thus allows us to clearly see some features of the evolution which are 
much more difficult in real QCD. 

We start in Section 2 by a short review of the elements we would like to discuss. 



2. JIMWLK and other letters 



We start by describing the general setting of the high energy evolution. We consider the leftmoving 
projectile with the wave function \P) colliding with the rightmoving target with wave function |T). 
The total rapidity of the process is Y, of which the target carries rapidity Yq. The S-matrix for 
the process is given by 



S(Y) = j Da a T [r] W^[a T (x,r)) ^_ Yo [a T (x, r)] . 



(2.1) 



Here a denotes the gA~ component of the vector potential due to the target charge density pp. 
It defines the scattering matrix of a single gluon projectile at transverse position x on the target, 
which in the eikonal approximation is given by 1 



S(x;0,t) = V exp{i / dr'T^^x, r')} . 
Jo 



(2.2) 



where T£ c = if abc is the generator of the SU(N) group in the adjoint representation. The field ay 
obeys the classical equation of motion || [jj 

a T (x,r)T a =g 2 ±(x - y) {st(y;0,r) p T a (y,r)T a S(y;0,r)} (2.3) 



The function S p in eq.(2.1) is the eikonal scattering matrix for a composite projectile averaged 
over the projectile wave function, see [Efl 



S p [a T ] = J dp P W p [pp] ex p{^ dT J d 2 xp a P (x,T)a^(x,T)^ (2.4) 



with Xi - the transverse coordinates of the gluons in the projectile. W p is the "weight functional" 
which determines the distribution of the charge density in the projectile. 

The evolution of the S-matrix with energy is given by the general expression of the form: 

d f $ 

— S = I Daft W$ n [a T (x,T)] X a T ,- 



S = f Da a T W? [a T (x,T)} X 



5 ax 



T l p -_ Y [a T (x,T)] . 



(2.5) 



The kernel of the evolution \ can be viewed as acting either to the right on S or to the left on W: 
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W T [a T ] ■ (2.6) 



In our convention the variable r is rescaled to run from to 1. Thus r = corresponds to i = — oo while 

T — 1 to X~ — +00. 



Note that according to eq.( |2.4D , the projectile- averaged S-matrix E p is the functional Fourier 
transform of the projectile weight functional W p . Thus ax plays the role of conjugate to pp. The 
evolution of S p determines the evolution of its Fourier image W F : 



by 
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(2.7) 



Since the functional form of the evolution of the target and projectile weight functionals should be 



the same (projectile-target democracy - PTD) comparing eq. (|2.7D with eq. fl2.6|) we conclude that 
the kernel x must be selfdual 
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(2.8 





5 ' 




\ ' 


X 


a, — 
da 


= X 


-1>T-,IP 
. S P _ 



where now a and p both refer to the same object, either target or projectile. Thus if the target 
and wave function are treated on equal footing, the kernel x must be invariant under the Dense 
Dilute Duality (DDD) transformation pf 



p(x, x 



5a(x, x 



5p(x, x 



ia(x,x ) 



(2.9) 



So far the approximations used in the literature to derive the kernel do not treat the target 
and the projectile equally. The approach |], f?j assumes that the the color charge density in the 
target is parametrically large (p a = 0{l/a s )). This yields the the so called JIMWLK evolution 
operator || [?J 
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(2.10) 



5a a (x,0) 5a a (y,0) 5a a (x,l) 5a a (y,l) 



5a a (x,0) 



S(z;0,l 



5a b (y,l) 



In this expression the functional derivative in eq. ( |2.10| ) do not act on S in the kernel, but on 
external factors only. The same result follows from the calculation || which assumes that the 
projectile carries only small color charge density. The derivation of |^4j on the other hand assumes 
that the target is dilute (small charge density). This calculation yields 



X 



KLWMIJ 



(z - x)i(z - y)i 



271-2 Jx,y,z (Z ~ X) 2 (Z - y) 2 

{p a (x, 0)p a (y, 0) + p a (x, 1) p a (y, 1) - 2 p a (x, 0)R(z, 0, l) ab p\y, 1)} 
where the "dual Wilson line" R is defined as 

r — n a b 

t f x rlz-T c - 6 

R(z;0,x~) ab - ^ Jo 1 



(2.11) 
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(2.12) 



Just like in eq.(|2.10|) the functional derivatives in R act only on external factors and not on the 



explicit factors p a (0) and p a (l) in Eq. (2.11). Neither one of these expressions is self dual, rather 
they are related by the DDD transformation eq.(|2.9|). 



Note that R acts as the finite shift operator. When acting on a function of p it shifts the 
charge density by the charge density of a single gluon 



R(z,0,l)^[p a (x,x-)] = Z y [p a (x,x-)+T a 5 2 (x-z)] . (2.13) 

More recently a step forward has been made by resumming part of the corrections due to finite 
density effects in the wave function of a more dilute object (projectile in the case of JIMWLK, 



and target in the case of KLWMIJ) [p5|, 26, 31, The reulting improved expressions are 



JIMWLK+ 

A, 



2 6^(z,0,[A])^ (z ,o,l)^(z,l,[A])| (2.14) 



KLWMIJ+ 



■^Jm*A\p])bU*A\p]) + b a (z,o,[ P ])b a (z,o,[p}) 

2bUz,0,[p])R ab (z,0,l)bUz,lM)}- (2-15) 



The field bf satisfies the "classical" equation of motion 

{di + 9 f ab %(z, [p])} d+b^z, [p]) = gp a (z, X-) ■ 

e^dib^z, [p]) - d 3 b*(z, [p]) + gf ahc b\{z, \p])bftz, [p])} = . (2.16) 
Expressions eqs,( [2.14| ) and ( p. 15 ) are again related by the DDD transformation eq.( |2.9| ). 



The KLWMIJ evolution equation was derived in |24|] by directly considering the evolution of 
the projectile wave function. In this framework the charge density operator depends only on the 
transverse coordinates p a {x). As quantum operators the charge densities at the same transverse 
coordinate do not commute 

[p a (x),p b (y)} = if abc p c (x)5 2 (x - y) (2.17) 

Thus in general the averages of powers of p can not be represented simply by a functional in- 
tegral with a weight functional which depends on the classical field p a (x) only. This is unlike 
JIMWLK, where such a representation is possible. It is however possible to have a functional 
integral representation if the " classical field" is endowed with an additional " ordering" coordinate 

{p a H Xl )...p a "(x n )) = J dp(x,x-)p a '(x 1 ,x^)...p a "(x n ,x-)W T [p(x,x-)] (2.18) 

On the right hand side of this expression the values of the additional ordering coordinate inherit 
the ordering of the operator factors on the left hand side: x^ > x^ > ... > x~. Clearly the weight 
functional W can not be arbitrary. It must ensure two basic properties of the " classical" correlation 
functions. First, the correlation function can not depend on the values of the coordinates x~ , but 



only on their relative ordering. Thus the correlation function should be a piecewise constant 
function of x~ . Second, under the interchange of the two adjacent values of x~ the correlation 
function must emulate the effect of the commutator eq.( |2.17 ) 



D[p}W[p}...p a (xr)p\ x - 



D[p]W[p]...p a (xj)p b (xr)... + if 



eabc 



D[p]W[p]...p c {xT)... 

(2.19) 

In particular this last property means that W can not be real, and thus cannot be given the 



meaning of probability density. It was pointed out in [24] that a similar problem has appeared 
in the literature before [38| where it was shown that these properties are satisfied by a weight 
functional which contains an exponential of a Wess-Zumino term. 



3. On the operator ordering, x vs x + and the Wess-Zumino-less weight func- 
tional. 

In this section we make several remarks on the structure of the evolution described in Section 2. 
3.1 On intimate relation between target and projectile 

A somewhat curious feature of the JIMWLK equation which is not frequently stressed in the 
literature, is that the same equation arises as the result of two apparently unrelated derivations. 
The original derivation H was given directly for the target wave function (equation for W T 
in eq. ( |2.6[) ), assuming explicitly that the target is dense, namely that the average of the single 
gluon S matrix is close to zero. Nothing has been explicitly assumed there about the properties 
of the projectile. On the other hand, one can deduce the kernel ^ JIMWLK i n a different way 



assuming the projectile is small and deriving the evolution of the projectile wave function [36, 25]. 
This then yields the evolution equation for the target Wt via the functional Fourier transform, cf. 
eq.(eq. ( |2.6| )). In this derivation we apparently do not assume anything about the properties of 
the target, instead we have assumed that the projectile is dilute. Nevertheless the two ways lead 
to the same evolution equation for Wt- The question is how come one obtains the same equations 
under seemingly unrelated conditions? 

The answer is that these two conditions - diluteness of the projectile and dense nature of the 
target are actually inextricably related through the form of the 5-matrix eq. (|2.4| ). The JIMWLK 
kernel gives a correct approximation for the full kernel of the evolution of Wt in the limit when 
«T is large, or equivalently when S is close to zero. On the other hand the evolution of S eq.( |2.5| ) 
is given by the functional integral over all ax- For W T which is sharply peaked around S = 0, the 
integral picks out the large ay region of S^fay]. Thus the approximation is safest when S p [ar] 
is also peaked around large ay. 

This however is precisely the condition for the projectile to be dilute. This is simply seen 
from eq. (2.4). Mathematically it is the general nature of Fourier transformation that E p [ot] is 



peaked at large ax if its Fourier image W p is dominated by small values of pp. The diluteness of 
the projectile is simply the statement that the projectile wave function contains only Fock space 



components with small number of gluons. Parametrically, if ay is of order one as is assumed in 
the derivation of H the projectile charge density that is picked out by the Fourier transform 
is of order O(l) also, which corresponds in our normalization to 0(1) number of particles in the 
projectile. In this case the projectile averaged S"-matrix T, p has only several first terms in the 
Taylor expansion in powers of S, as each projectile gluon contributes a factor of S. Thus for a 
dilute projectile, the Taylor series expansion of T, p has only small number of terms and is indeed 
dominated by S ~ 0. 

The diluteness of projectile, and the dense nature of the target therefore go hand in hand. The 
corollary to this observation is that if one wishes to discuss scattering in the situation where the 
projectile is not dilute, the JIMWLK approximation to the target weight function is not necessarily 
applicable even if the target is dense. The projectile averaged S matrix T, p is then itself peaked 
around S ~ 1, and so even if the target is mostly dense, the tail of the distribution in W T at 
ax C 1/g will contribute significantly to the total cross section. Physically this is very reasonable. 
If the projectile is dense, it has large probability of scattering on a variety of configurations in 
the target wave function: the ones with large charge density as well as the ones with low charge 
density. Thus the tail of the distribution Wt which determines the weight of dilute configurations 
in the target wave function becomes important. A discussion of the evolution of the total scattering 
amplitude in this situation is to some extent superfluous, as it vanishes at all energies. However 
there are other less inclusive quantities which we believe will be strongly affected by this tail of 
the distribution. An example of such a quantity is diffractive scattering which is sensitive to the 
dilute peripheral region of the target. The single gluon emission spectrum also may get sizeable 
contributions from the target configurations in the tail of the distribution. 

3.2 x + or x~1 

Our second remark is essentially semantic. However since there has been some lack of clarity on 
this point in the recent literature, we thought it worth a comment. 

The issue in question is whether to call the ordering variable entering the KLWMIJ kernel x~ 
or x + . We prefer simply to think of it as of an ordering variable and clearly in practical terms its 
name does not matter. There is however a certain physics picture that suggests its analogy with 
x~ . In particular consider the KLWMIJ+ equation eq.( p,15 ). It involves the factors depending 



on the "classical" field 6j(x,r), which satisfies an equation of motion identical to the classical 
Yang-Mills equation in the presence of the source p(x,r), where t is the ordering variable. In 
our wave function approach [^4|] this equation arises in the following way. To begin with one 
has to solve operator equations for the operator bi(x) which has an operator source p(x) without 
any r-dependence. Solving these equations, for example in the perturbation theory, one has to 
keep track of the ordering of operators p a as explained in the previous section. This procedure is 
streamlined by introducing the ordering variable r and treating all fields as classical. The resulting 
" classical" equation is identical to the classical equation solved in |?7]] if one identifies the ordering 
coordinate t with the longitudinal spatial coordinate x~ . It is actually more consistent to interpret 
the ordering variable as the rapidity lnx~/xQ rather than the longitudinal coordinate itself, since 
it is dimensionless, but we will keep referring to it x~ for simplicity. This identification is natural 



in the sense that x~ and not x + is the spatial coordinate of the target. As the whole action of 
KLWMIJ takes place in the projectile Schroedinger (light cone) wave function there is no room 
for the time coordinate in this framework. 

Once we have identified the ordering coordinate of 6, as x~ , we have no choice but to identify 
the ordering coordinate of R in eq. (|2.12j ) well. This is quite clear, the operator R simply 

shifts the charge density of the target by the charge of an extra gluon ( |2.13| ). Again, since this is a 
gluon in the light cone wave function, its coordinate is naturally identifiable with the longitudinal 
coordinate rather than the light cone time. 

In fact this interpretation goes hand in hand with the usual spatial picture of the low x gluon 
field. One normally thinks of the valence gluons as concentrated in a "pancake" of small extent 
in x~ , while the wee gluons are distributed in a wider longitudinal slab at smaller values of the 
spatial longitudinal coordinate z, and so at smaller values x~ . Eq.( 2.15 ) conforms nicely to this 
picture. In eq.( |2.15| ) the values of x~ in the various factors are ordered from left to write. Thus the 
values of x~ in the factors of the charge density in the leftmost classical field bi are the smallest, 
the values of 5/5p(x~) in the factor R are larger and the values of x~ in //s of the rightmost bi are 
the largest. Following the derivation of 24] we can think of the rightmost bi and R as belonging to 
the wave function, while the leftmost bi belonging to the conjugate wave function. The rightmost 
bi in eq.( 2.15|) corresponds to the valence gluons in the wave function and thus takes up larger 
values of x~ . The wee gluons are the ones created by the shift operators R, and consistently with 
our interpretation take up smaller values of the longitudinal coordinate. From the action of R 
eq.( |2.13| ) it is obvious that the gluon emission probability is independent on x~, which is again 
consistent with the common intuition that the distribution of the wee gluons is flat in rapidity. 

The weight function W on which the action of the kernel is defined therefore depends on the 
field p(x~), where x~ takes values in three intervals necessary for representation of both factors 
of bi and R. This contour is visualised in Fig. |l]. The same structure is true for the JIMWLK+ 
equation if it is understood as the equation for the evolution of the target wave function in the 

1 . Of course if one writes 



regime of dense projectile. This is the logic adopted in Refs. [24, 26 



the equation for the evolution of the projectile weight function the natural identification of the 
ordering coordinate there would be x + . 
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Figure 1: Ordering structure of (a) JIMWLK+ kernel, (b) KLWMIJ+ kernel, and (c) W 



The expression that is at the first sight somewhat puzzling within this interpretation, is 
eq.(|2.4D, as there both ax and pp carry the same ordering variable. However one should re- 



member how this equation arises [26]. The operator expression for the eikonal scattering matrix of 



the projectile charge density of the target field is 

t p =V ex p|^ dx ~ J d 2 xp a P (x)a^(x,x-)^ (3.1) 

where the path ordering is along x~ . The natural variable here is x~ , since it is the spatial 
variable from the point of view of the target fields and time coordinate as far as the projectile is 
concerned. At this point pp is an operator and does not carry any longitudinal label. To average 
this expression over the projectile wave function we should expand eq.( |3.l|) in powers of pp, endow 
Pp with an ordering coordinate, which we would now call x + (as it is the spatial variable in the 
projectile wave function), and average with the weight function W[pp(a; + )]. We observe however 
that since eq.( |3.1| ) is path ordered with respect to x~ , the ordering of the operators pp in any term 
in the expansion is the same as the ordering of x~ coordinates in the factors olt which multiply 
them. Thus the values of variable x + are always ordered in exactly the same way as the vaues of 
x~. Since the correlators of pp(xf) depend only on the ordering of x^'s, we can simply set the 
value of the longitudinal coordinate of pp equal to that of x~ in olt- If one wishes to find a physical 
picture behind this mathematical statement, the proper words are that the projectile charges and 
the target fields always interact at the point t = 0, where x + = x~ . Thus at the interaction point 
indeed the longitudinal labels of the projectile density and the target field are equal. 

3.3 Life without Wess-Zumino. 

Though the explicit form of the WZ term is known (see, e.g. |24|| ), it is not clear if this knowledge 
can be of any use to perform the actual averages for physical observables. In this remark we show 



that the following compact parameterization of the weight-functional W [28] 

W[p(x,x-)} = E[R]5[p(x,x-)} (3.2) 
with £ - an arbitrary functional of R, does not need to be augmented by the WZ term. This 



form is a generalization of the form used in [13] in the framework of the dipole model with the 



difference that it explicitly incorporates the ordering variable, and the expression for R should not 
be truncated at low order in derivatives as in [|l3| ]. 

The logic for this "ansatz" is the following. Consider a hadronic state containing energetic 
gluons in its wave function. When scattering on any hadronic target, its S-matrix must be some 
function of the eikonal factor for the propagation. As the S'-matrix is given via eq. (|2.4j ) we must 
have 

/ Dp a e l f pa(x ' T)aa{x ' T) W[p\ = Z[S] (3.3) 

Eq.(|3.2j) is specifically devised to satisfy this requirement and thus is a natural representation for 
a weight functional. We will now show that it also automatically satisfies the requirement eq. ( |1 . l\ j 
and thus does not have to be augmented by a Wess-Zumino term. The x~ ordering is relevant for 



p's which are at the same transverse coordinate, thus we omit the transvers coordinate dependence 
in the following. 



W will use the fact that p acts on R as rotation [28]: 

P a (x-)R bc (yI,y2) = [R{yux~)T a R{x- ,y^)] bc efaj - x~)@{ X - - yr) (3-4) 
To illustrate the general pattern, we consider first the one and two point functions 



(p a {x-)) 



J d PP a { x -) X[R]8[p] = J dp [R(0, 



x -)T a R(x-,l)\^^^8[p] 



Jcd SR dc 



,a SE[R] 



- cd 



SR, 



\R=1 



dc 



(3.5) 



(p a (x^)p b (x 2 )) = I dpp a {x^)p b {x 2 ) E[R]S[p] 
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(3.6) 



The correlator is clearly piecewise constant, and depends only on the ordering of the coordinates. 
Moreover, the last term is symmetric under exchange of x\ and X2, while the interechange in the 
first term yields precisely the one point function, as required by eq.( |l.l| ) 



(p a (x^p b (x^)) = (p a (x^p b (x^)) +;r fec <p c (* r )>sign(:r 2 - - x^) 



(3.7) 



Quite obviously the same holds also for all higher correlators. The result for the three point 
function for example is 

(p ai (^)p a2 K)p a3 K)) = / dpp(x^p(x^)p c (xz) E[R}8[p] 
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+ [T a 1T a 3T a 1]cdQ{x - _ X -)Q( X ~ _ x ~) + [ T a ^ T a % d Q{ X ^ - X 3 ) @(x 2 - X^) 

+ [T a * T a * T a % d Q{ X2 - - X3)0(xi - x 2 ) + [T a2 T a3 T ai ] cd e(xJ - xj)0(x]f - x^)} 



1 r» 



5E[R] 



5R. 



dc 



R= i + {[r ai r a2 ] cd rf;e(x- - x~) + [T a * T a % d T^ } ®{ x - 



+ C [r-r a 3] e/ e(x 3 - - x 2 -) + r- [rr] e/ e(, 2 - * 3 -)} ^§H«=i 

Again we see explicitly that eq.(l.l) is satisfied. This illustrates the general pattern: the factors 
p ai (x~ ) in the correlator are replaced by the color matrices T ai ordered according to the ordering 
of x~ . That is for the correlator (p ai (x^) . . . p a% {x~) p a "(x~)} with strictly order 
... < x~ we get 



n 



(p ai . . . p an ) = ^ Yl t Tai • • • T a ^] cd [T a (n+D . . . T a '2] e/ • • • [T a ^-i +1 > . . . T an ] bg 

k=l l<u<«2<---<ifc-i<n 



5R dc ■ ■ ■ 5R gb 



\R=1 



(3.9) 



For different orderings of Xi , the n-point function is obtained from eq.( |3.9D by the rule eq.(l.l 



In the dense regime the introduction of the longitudinal coordinate is not necessary || |2^] 
and the correlators must become symmetric under the interchange the factors of p. This should be 
seen directly from eq.( |3.9| ). Indeed the wave function of a dense system is characterized by a large 
number of gluons, say K ^> 1. In this case S[i2] behaves roughly as R . The n-point function (as 
long as n < K) is dominated by the n-th derivative of E, which as we have seen above, does not 
depend on x~ at all: 

{P ■■■ P ^ 1 c d ---h 9 SRdc . . . 5R gb lR=1 (3 - 10) 



In accordance with discussion in the previous sunsection, we note that the variable R in eq.( |3.2|) 
is defined as the path ordered product along all three segments of the contour in Fig.|l]. That is we 
must understand R in eq. ( |3.2| ) as the product R = RiRjiRju (Fig. |l|,c). This does not affect the 
derivation of this section, and also correctly reproduces the correlation functions involving factors 
of p that appear due to expansion of the fields bi[p] in KLWMIJ+, eq.(2.15). 



4. Self-duality and the dipole toy model. 

A simplified version of the high energy evolution - the dipole model was introduced by Al Mueller 
Q. It describes the leading high energy behaviour in the large iV c limit as long as the densities 
in the wave functions are not too large. As shown in [^] the dipole evolution equation can be 
obtained as a well defined limit of the JIMWLK evolution. We can define the projectile dipole 
creation operator s and the target dipole creation operator r: 



s(x,y) = Tr[S' F (x) S F (y)] ; 



r{x,y) = Tr[W F {x)R F {y)] 



(4.1) 



where F indicates the fundamental representation. The canonical annihilation operators are simply 
and 4- for projectile and target respectively. 

If the target weight function is a function of s only, that is W T = W T \s\, the action of the 
JIMWLK kernel on it in the large N c limit is equivalent to the action of the dipole kernel |H| (see 
also p|]): 



X 



JIMWLK 



W T [s] 



dipole 
Xdp 



Ss 



W T [s] 



(4.2) 



The dipole kernel in the dilute projectile regime can be obtained by reformulating the original 

5 



Muller's model and it has the following form [15] 

5 



dipole 
^-dp 



Ss 



K 



s(x,y) - s(x,z)s(y,z)] 



x,y,z 



5s(x,y) 



(4.3) 



On the other hand for dilute target which contains only dipoles in its wave function W 1 



W T [r], and the action of the KLWMIJ kernel in the large iV c limit is equivalent topi]: 



X 



KLWMIJ titT 



W 1 [r] = x7t 



dipole 



' Sr 



W T [r] 



(4.4) 



Just like the JIMWLK and KLWMIJ kernels, the dipole kernels in the two regimes is related 
by the duality transformation eq.( |2.9D . Since the dipole model is much simpler than the complete 
evolution equation, one may hope to gain some insight into some features of the "holy grail" self 
dual evolution by studying it. So far unfortunately there is no exact formulation of duality in the 
framework of the dipole model. We are unable to provide such a formulation here, instead we will 
consider an even further simplified toy model of the evolution, which does not involve transverse 
coordinates. Physically this corresponds to considering the evolution of distribution of dipoles of 
fixed size at a given impact paramter. 

The projectile averaged S matrix S p of the toy model is given by 



(4.5) 



Here s is the scattering matrix for the scattering of a single dipole on the target and P n is the 
probability to find n dipoles in the projectile. This probability plays the role of the projectile 
weight function 

W p (n) = P n . (4.6) 
In the dilute projectile regime the evolution of S p is given by 



dY 



with the dipole kernel X d % ^° le 



dipole 



ds 



dipole ^ I 
Xdp ^ 



ds 



(4.7) 



(4.8) 
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where we can think of ^ as of the dipole annihilation operator, while s as of creation operator. 
The evolution ( |4.8| ) is equivalent to the evolution of the probabilities P n [15] 



dP d 

= a s (n - l)P„_i - a s nP n = a s (e~*r - l)nP n (4.9) 

This is a simple hierarchy of equations which has an interpretation that in every step in the 
evolution one additional dipole is emitted with probability a s . Note that n is a discrete variable 
and thus the action of d/dn is not well defined. However the action of the shift operator e~*r on 
any function F n is perfectly well defined and gives 

e ~£ F n = F n _i (4.10) 

Thus the analog of the dipole version of the KLWMIJ equation is 

Xdt P ° le ( n > e ~ 3 ") = ®s(e~~^ - l)n (4.11) 

The shift operator e~*r is the analog of the "dual Wilson line" R and creates an extra projectile 
dipole. 

Following Ref. 0] we can define the dipole densities p\: 

The evolution eq. ( |4.9| ) or its equivalent eq. ( |4.3| ) leads to the equation for pP k |2(| 

j± =a s {kfF k + (k- 1)^_J (4.13) 

It is sometimes useful to define the amplitude of simultaneous interaction of k dipoles on the 
target - 7^. The S'-matrix can be written in terms of pk and jk as [|| 

S(y) = 1 + ^(-l) fc ^ 7fc (4-14) 

k 



As shown in [20], the evolution eq. (|4.13p can be mapped onto the Balitsky's hierarchy for the 
amplitudes 7: 

^ = a s k [ 7fe - 7 fc+ i] (4.15) 

It is instructive to see how the duality transformation works in this simple model. It is useful 
to keep in mind the following set of anologies (for the definition of a see below): 

p -> n; a^nlna; S^e nlna ; R -> e~ & 

The S-matrix is computed by averaging T, p over target distribution Ty T [s] 

S(Y) = j ds^_ Yo [ S ]Wl[s] (4.16) 
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Written in this form W T gives a distribution of probabilities that a single projectile dipole scatters 
off the target with the amplitude s. Though we have written integral over s, in fact s is discrete, 
and the integral over s simply should be understood as the sum over postive integers n with s = a n . 

Let us denote by a a cross section of a single projectile dipole to scatter on a single dipole of 
the target (an analog of two gluon exchange in real QCD), a = 1 — 70^ with 7 of order unity. If 
we allow a projectile dipole to multiply rescatter on m dipoles of the target, the natural model 
for the scattering amplitude on a m-dipole target state is s m = a m . The target by itself has a 
distribution of dipoles 

W T (m) = T m (4.17) 
Then eq. (4.16) can be written in a manifestly symmetric form 



S(Y) = ^P n (Y -Y Q )T m (Y )a mn = £ P n T m e mnl ™ (4.18) 

n,m n,m 

To derive the duality transformation in our toy model we assume that P n evolves according 
to the action of some kernel x dipole ( n , e~*r) acting on P n : 

§pi = X dip ° le (n, e"*) P n (4.19) 

By "integrating by parts" we take the kernel x p ' to act on a nm . The result of this integration 
by parts is that d/dn — > —d/dn, and the ordering of factors containing n and d/dn is reversed. We 
then notice that 

ne^a nm = e m lna — ^_ a nm ( 420 ) 

In a dm 

Integrating by parts again we take the kernel to act now on T m . Th.6 result of this three step 
procedure is the evolution equation for T m 

i>[ "' X diP ° le (-T^-T-^ mhla ) T m (4.21) 



dY \ ln<7 dm 

where the function \ differs from the function x by the opposite ordering of the factors containing 
n and d/dn. In order to obtain eq. (4.21) we used the Lorenz invariance (LI), which requires that 



iS does not depend on Yq. Now the projectile-target democracy says that P and T have to obey 
the very same evolution law. Hence we deduce 

dipole ( ' e -4-\ _ -dipole ( \ ^ e « lnc 



dtpole 1 e -^\ = ^pole " „ in CT u ^ 

V / y ma dn ) 

The DDD transformation can therefore be formulated as the transformation 

& -k -> e nln ° (4.23) 

In a dn 

accomapnied by the interchange of the order of n- and <i/<in-dependent factors. 

We remark that the DDD transformation ( 4.23| ) as well as the selfduality condition ( |4.22| ) have 



been derived in a specific model for the projectile-target interaction, namely (e mn ln<J ). Although 
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this is a reasonable ansatz, in our toy model we are not bound to choose this particular form. Have 
we choosen another ansatz for s the resulting transformation would be different. Thus although the 
existence of duality transformation follows from the Lorentz invariance and PTD, its mathematical 
action depends on the form of the projectile-target interaction. 

Let us now apply the DDD transformation ( [4.23| ) to the "dilute" kernel eqs. (4.8,4.9). The 
resulting kernel is analogue of the JIMWLK equation. 

To make sense of the operator -j- we should represent it in terms of shift operators exp{— k d/dn}. 
This reresentation is obtained using the trivial property x = — lne~ x : 

k=l 

The evolution equation then can be written entirely in terms of the creation operators of k dipoles 
e~ k 3n (this interpretation is in analogy with eq.(f4.9j)). 



dipole 

Xdense ~ a s A 

k 



The expansion ( [4.25 ) has a remarkable propery: each term in the expansion separatly conserves 



probability. The subtraction of the unity takes into account virtual contributions for /v-dipole 
emission 2 . 

The eq. ( |4.26j ) suggests the interpretation that in the dense regime, arbitrarily many dipoles 
are created during one step of evolution. This interpretation is indeed correct, and is consistent 



with the actual evolution in QCD [39, 12]. However it is also clear that the kernel does not correctly 
account for the propagation of these newly born gluons through the target, as only one eikonal 
factor is present in eq.( [4.24 ). This is an analog of the situation discussed in the beginning of the 



previous section. Assumption of the dense wave function of one of the colliding objects (multiple 
powers of exp{— 4tY) goes together with the assumption of the dilute wave function of the other 
(only first power of exp{ n lncr}). In the full kernel higher powers of the dipole creation operator 
must be also accompanied by higher powers of the eikonal scattering amplitude. The following 
interpretation of eq. (4.26) looks natural. Real emission of k dipoles in the full evolution kernel 



is accompanied by virtual corrections to all wave function components with numbers of dipoles 
smaller than k. The component of the wave function with k emitted dipoles contributes to the 
S-matrix with the factor s k = a nk , while the virtual terms contribute with lower powers of s, 
corresponding to the number of dipoles in the state whose probability is being "repaired" by the 
virtual correction. In the dense limit s < 1 and so the leading contribution to the S'-matrix comes 
from the propagation of the state with only one emitted dipole and the state with no emitted 
dipoles. Thus although only one emitted dipole shows up as the "real" emission term, all the 



2 One should be able to obtain a similar representation for the JIMWLK equation by representing J- as an infinite 



series in powers of (R— 1). 
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higher dipole contributions show up in the virtual correction piece: the norms of the "no extra 
dipole state" and "one extra dipole state" are decreased by the total probability of emitting more 
than one dipole in one step of the evolution. This interpretation is further strengthened by the 



observation that the series in the one dipole creation operator in eq. (4.26) is an alternating sign 
one, which is what one naturally expects from a virtual correction. 

In the dense limit where n is large and ^ is small we can actually simplify the kernel by 
keeping only the first, single dipole term in the expansion ( 4.25j ) . 

x dipole = - g J_ ^ _ ^ _ e « Ina^j (427) 

The evolution eq. ( |4.27 ) possesses a remarkable property: it is selfdual under the DDD transforma- 
tion ( 4.23 ). The kernel eq.( 4.27| ) is expressed entirely in terms of the "Wilson lines" and in the two 
known limits of either dilute target or dilute projectile reduces to the "JIMWLK" and "KLWMIJ" 
equations. The kernel eq. ( [4.27 ) leads to the following evolution equation for the probabilities 



dY 



lna\ 



[(1 



a 



-v. 



n-1 



{l-T n )Pn 



(4.28) 



This has an interesting interpretation. The coefficient in front of P n on the right hand side has the 
meaning of the probability per unit rapidity to emit an extra dipole from the n-dipole state. In this 
respect eq. ( |4,9| ) is clearly problematic, since for large n the emission probability grows without 
bound. Therefore it can only be considered informative for relatively small n, such that a s n < 1. 
The behavior of eq. ( [4,28| ) is much better, since the emission probabilities are finite even for large 
n and so the gluon emission probability "saturates". Since by definition a < 1, the limiting value 
of the probability in eq. ( 4.28] ) is a s /|lna|. We note however that this limiting value generically 
is still larger than unity. Within our present model a s is in principle unrelated to a. However 
in QCD the analog of a is the dipole-dipole scattering amplitude, which behaves as a = 1 — 70^ 
with 7 of order unity. With this parametric dependence the limiting value of the coefficient in 
eq. ( 4.28j ) is 0(l/a s ) ^> 1. Thus although the emission probability given by the selfdual form of 
the evolution saturates, it still has the probabilistic interpretation only for a s n < 1 just like the 
initial expression eq.(|4.9|). 

The expression eq. ( 4.2§| ) is analogous to the self dual variant of JIMWLK+ evolution recently 
proposed in Ref. [|31jJ. The lesson we learn from our toy model is that self duality is not sufficient 
to gurantee the accuracy beyond the asymptotic limits large and small density limits. In general 



we should expect higher powers of both R and S to enter the fray (see also Ref. |34|] for discussion 
of the selfdual form of the QCD evolution). 



A somewhat surprising aspect of the evolution eq. (|4,28|) is that just like eq, (|4.9| ) it describes 
only "splitting" of dipoles and does not contain "mergings" in the parlance of [22, 21, |23]| . One 



could naively expect to find merging terms which lead to growth of the probability P n proportinal 
to probabilities P^ with k > n. Such terms do not appear in eq.( 4.25| ) nor in the self dual limit of 



the evolution eq.(4.27). On the other hand the common wisdom is that the "splitting" evolution of 
the projectile dipoles is somehow equivalent to "merging" evolution in the target. The resolution 
of this paradox is that even though there are no merging terms in the evolution equation for the 
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probabilities P n , they do indeed apear in the equation for the n-dipole densities p n and the n-dipole 
scattering amplitudes 7 n . In particular if instead of canonical dipole model eq. (|4.9|) we start from 



the selfdual evolution eq. (4.28) we obtain the following evolution in terms of pj. 



dY 



In cr 



v 

Pk- 



£ fP n £ (In arC(k,m,n) 



n=k— 1 



m=n+l— fc 



(4.29) 



where the coefficients C(/c, m, n) are computable combinatorial factors. Since the expansion in lncr 
is essentially an expansion in truncating it at second order we obtain 3 



dY 



O.s- 



k(? k [l + (k--) In a] + — 1) Ph-i [1 + (^ 



. . k(k + l) . p 
lna] + - lnapj^ 



(4.30) 

The first two terms are the same as in Balitsky's hierarchy amended by oq corrections. The last 
term corresponds to the merging of Refs. 22, 23|. When translated to the equation for the 
amplitudes 7&, this term leads to a contribution of the form 



d-jk 
dY 



~ a. 



|lna|7fc_i 



(4.31) 



which is generally understood as Pomeron merging. Thus curiously enough we see that "Pomeron 
merging" term in the evolution of the amplitudes 7^, when translated into the language of proba- 
bilities P n is in fact no different from dipole splittings. Physically this is perfectly reasonable, as 
boosting the hadronic state can never lead to the decrease in number of gluons or dipoles in the 
wave function. These numbers always grow, but the rate of their growth decreases at high energy. 
The Pomeron mergings describe the decrease in the rate of growth of the cross section, and not 
the decrease in the number of partons. 
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